METHOD AND APPARATUS FOR DETERMINING PARTICLE SIZE AND 

COMPOSITION OF MIXTURES 

This application claims, priority under 35 U.S.C. §119 (e) 
from provisional patent application serial no. 
60/449,418 filed on February 24, 2003, the entire 
content of which is incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to the quantitative 
determination of analytical properties of chemical 
5 systems . More particularly to a method and apparatus for 
using broadband, preferably near-infrared light to 
determine both composition and particle size of colloidal 
mixtures. The present invention also relates to a method 
and apparatus for using narrow-band or laser light with 
10 autocorrelation techniques for characterizing variations 
in the properties of creams and lotions. 

2. Prior Art 

The classical method for obtaining absorbance spectra is 
via direct light transmission through a transparent 

15 saxDple. A polychromatic beam suffers a loss of intensity 
at particular wavelengths as it traverses a distance L, 
known as the path length, through a sample. Absorption 
of light is a consequence of interaction with the 
molecules of the sample, particular species favoring 

20 specific wavelengths of light. By calculating the ratio 
of the transmitted beam to the initial incident beam, one 



obtains the transmission spectnua. Taking the negative 
logio of the transxoission spectrum provides the more well- 
known absorption spectrum. In general, these systems are 
well described by Beer's Law, which shows that the 
5 absorbance is proportional to the path length. 

Diffuse reflectance spectroscopy is the extension of this 
method to non-dilute, turbid, and even opaque samples. 
For these types of saaxples, light is backscattered, i.e., 

10 directed back in the same direction as it entered. Soblb 
of the light may have been absorbed as it penetrated the 
sample, analogous to phenomenon in transmission. In 
these systems, however, the path length is determined by 
the amount of scattering the light has undergone. The 

15 usual theory governing such systems is that of Kubelka- 
Munk. 

Recent work applying diffuse reflectance to physiological 
systems has made use of modified measurement techniques 

20 while expanding the utility of these measurements. For 
example, the work of Schmitt and Kumar, Spectral 
Distortions in Near- Infrared Spectroscopy of Turbid 
Materials" Applied Spectroscopy, 50, 8, pp. 1066-1073 
(1996) , which is hereby incorporated by reference in its 

25 entirety, was based on a variable s^aration distance 
between the optical fiber that introduced light into the 
system and the one that collected it. The authors show 
that by departing from the traditional backscattering 
geometry of routine diffuse reflectance, they can detect 

30 variations of the absorption spectrum that might not have 
been expected. For example, they show that for a fiber 
separation of 6 mm, they could get an absorbance 



equivalent to a path length of 24 mm. They attribute 
this discrepancy to photon diffusion and demonstrate how 
to calculate the amount of absorption in such mixtures* 



5 The critical feature that Schmitt and Kumar demonstrate 
is the dependence of the photon diffusion, and therefore 
the distortion of the spectrum away from the diffuse 
reflectance result, on the size of the particles doing 
the scattering. 

10 Dynamic light scattering for particle size determination 
is usually done in dilute systems. Instead of a 
broadband source, laser light is used. The signal is 
autocorrelated to provide a measure of the fluctuation of 
the intensity as a function of time, usually on a 

15 microsecond - second time-scale. Particles undergo 
Brownian motion and the resulting dynamics impacts the 
intensity of the Rayleigh scattered light, generally 
causing a decreasing autocorrelation function which 
frecpiently appear to resemble exponential decay. When 

20 this signal is analyzed, diffusion coefficients, and 
inf erentially , particle diameters can be deduced. 

SUMMARY OF THE INVENTION 

It is an object of the invention to provide a method and 

25 apparatus for monitoring particle size in a process 

wherein particle size of a material varies as a function 
of time. 
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It is a further object of the invention to provide a 
method and apparatus that characterizes particle size by 
reflection of light from the particles. 

It is a further object of the invention to provide a 
S method and apparatus that characterizes particle size by 
varying more than one parameter in the apparatus 

The invention comprises a method of determining the state 
of a mixture comprising^ in general, both scattering and 
absorbing constituents, which conprises injecting light 

10 into the san^le with a fiber optic; collecting light from 
the sample with a second fiber optic; spacing the fiber 
optics at a variable distance from each other, the 
coupling of injection and collection fibers being called 
a probe, and coialyzing the reflected light with a 

IS spectroscopic analyzer. The variability of the spacing 
may be controlled by a mechanical stage driven by an 
electric motor. In turn, this motor may be operated by 
the computer that controls the spectrometer. 
Alternatively, the invention can be implemented by using 

20 multiple collection fibers strategically spaced with 
respect to the injector; using a multiplexing device to 
select the specific collection fiber to observe; and 
analyzing the reflected light from the selected fiber 
with a spectroscopic analyzer. 

25 

Such spectroscopic analyzer may be, but is not limited 
to, an FTIR, photodiode array, AOTF or grating 
spectrometer. Further, analysis involves comparing the 
suite of collected spectra with a comparable suite 
30 collected from a sample in the desired state (''finished 
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product") ; additionally or alternatively, calculating the 
average particle size or the particle size distribution 
based on the mathematical inversion of the photon 
diffusion equations; additionally or alternatively, 
5 calculating the average particle size or the particle 
size distribution from a database such as one that can be 
provided by autocorrelation methods (also known as 
^^dynamic light scattering") or other analogous methods 
(acoustics , etc • ) . 

10 

Highly concentrated, virtually opaque, systems such as 
creams and lotions do not meet the criteria for analysis 
by this method. The particles are moving too slowly 
since they are large and likely to be in a viscous 
15 matrix. The systems multiply scatter, further confusing 
the traditional interpretation of the scattering signal. 

Suspending the need to quantify a ^^partide size" has 
merit for characterization of such dense materials, 

20 however. In many cases a technician sixnply wishes to 
classify a material according to a series of known and 
reliable properties. Among these would be viscosity, 
color, sheen, and other variables known to have impact on 
the performance properties of the material. When the 

25 shape of the correlation function changes with sample we 
have a new means of making measur^ents in a class of 
materials heretofore unexamined by optical means . 

The invention is straightforward, noninvasive and 
economical . 

30 The invention is directed to a method for evaluating 

one or more materials in accordance with size of 
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particles therein, comprising evaluating a spectrum of 
light reflected from a first group of particles; 
evaluating a spectrum of light reflected from a second 
group of particle; comparing results of the evaluating 
5 of the first group with results of the evaluating of the 
second group; and providing an indication of a state of 
the material when the comparing produces a predetermined 
comparison result. 

The method may be used for monitoring a process wherein 

10 size of particles is changing as a function of time, and 
wherein the evaluating of a spectrum of light from a 
first group of particles comprises evaluating a spectrum 
of light reflected frani the particles at a first time; 
the evaluating of a spectrum of light from a second group 

15 of particles coniprises evaluating a spectrum of light 
reflected from the particle at a second time; the 
ccmqparing results comprises comparing results of the 
evaluating at the first time with results of the 
evaluating at a second time; further comprising providing 

20 an indication of a state of the process when the 
comparing produces predetermined comparison results. The 
evaluating may be done so that at at least one of the 
first time and the second time occurs during the process. 
The evaluating at at least one of the first time and the 

25 second time may be an evaluating of a reference sample at 
a time other than during the process. 

The method may further comprise varying distance between 
a first light conductor for conducting light to the 
particles and a second light conductor for conducting 
30 light from the particles for the evaluations, and 
evaluating spectra of light received from the particles 



at a plurality of distances between the first light 
conductor and the second light conductor. The distance 
may be varied to eliminate spurious reflections from the 
groups of particles or from a container in which the 
5 particles are disposed. 

The particles may be in a liquid, and selected from a 
group comprising at least one of solid particles and 
liquid droplets. 

The method may further cosqprise moving, at a plurality of 
10 different times, a first light conductor for conducting 
light to the particles and a second light conductor for 
conducting light from the particles along a column 
containing the liquid. The method may further comprise 
determining at least one of a rate of motion of an 
15 interface in the liquid, and characteristics of particles 
in a vicinity of the interface. The method may further 
comprise determining characteristics of particles on 
opposite sides of the interface. 

The method may also further comprising determining 
20 compositional characteristics of the particles by 
analyzing at least one of the spectra, the coxnpositional 
characteristics including at least one of particle size 
and chemical ccnaposition. 

The first group of particles and the second group of 
25 particles may be successive groups of particles in a 
batch of particles undergoing a process of mixing, and 
the evaluating of the particles may be performed at 
different times, to obtain an indication of homogeneity 
of a mixture containing the particles. 
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The first group of particles and the second group of 
particles may be successive groups of particles in a 
process flow stream. 

The first group of particles and the second group of 
5 particles may be successive groups of particles in a 
process, and the method may further comprise using the 
results of the evaluating to monitor changes in 
characteristics of the particles to obtain an indication 
of homogeneity of the particles. 

10 The method may be used to perform at least one of 
classification of materials, monitoring of a process, 
determining authenticity of a product, and determining 
quality of a product. The method may also be used to 
determine size of particle in at least one of the first 

15 group of particles and the second particles. The at 
least one group of particles may be suspended in a liquid 
characterized by a series of parameters and the size of 
the particles may be determined by utilizing a 
predetermined relationship between values of the 

20 parameters, a spectrum of reflected light from the 
particles and the particle sizes. The distance between a 
first light conductor for conducting light to the 
particles and a second light conductor for conducting 
light from the particles for the evaluations may be 

25 varied, and the utilizing repeated at a plurality of 
distances • 

The invention is also directed to aA method for 
determining characteristics of a material comprising 
illuminating the material with coherent light so as to 
30 produce scattered light; autocorrelating the scattered 



light; preprocessing a signal representative of the 
scattered light to produce a processed signal; performing 
at least one of a cximulant analysis and a cluster 
analysis to provide analyzed data; and associating the 
5 analyzed data with particular characteristics of the 
material . 

The invention is also directed to an apparatus for 
evaluating one or more materials in accordance with size 
of particles therein, comprising means for evaluating a 

10 spectrum of light reflected from a first group of 
particles and from a second group of particles; means for 
comparing results of the evaluating of the first group 
with results of the evaluating of the second group; and 
means for providing an indication of a state of the 

15 material when the comparing produces a predetermined 
comparison result. 

The apparatus may be used for monitoring a process 
wherein size of particles is changing as a function of 
time, and wherein the means for evaluating of a spectrum 

20 of light from a first group of particles coniprises 
apparatus for evaluating a spectrum of light reflected 
from the particles at a first tdLme and at a second time; 
the means for conparing results comprises comparing 
apparatus for comparing results of the evaluating at the 

25 first time with results of the evaluating at a second 
time; further comprising means for providing an 
indication of a state of the process when the conrparing 
produces predetermined comparison results. The means for 
evaluating evaluates at the first time and at the second 

30 time, during the process. 
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The ai^aratus may further comprise a storage device for 
storing a reference produced at a time other than during 
the process, for evaluating at at least one of the first 
time and the second time. 

5 The apparatus may further cCToprise means for changing 
distance between a first light conductor for conducting 
light to the particles and a second light conductor for 
conducting light from the particles for the evaluations, 
and meeuis for evaluating spectra of light received from 
10 the particles at a plurality of distances between the 
first light conductor and the second light conductor. 

The apparatus may further comprise means for containing a 
liquid in which the particles are dispersed and means for 
moving, at a plurality of different times, a first light 
15 conductor for conducting light to the particles euid a 
second light conductor for conducting light from the 
particles along the means for containing the liquid. 

The apparatus may further comprising means for 
determining a rate of motion of an interface in the 
20 emulsion, and means for determining characteristics of 
particles in a vicinity of the interface, as well as 
means for determining characteristics of particles on 
oj^osite sides of the interface or in a vicinity of the 
interface . 

25 The apparatus also may further comprise a flow cell 
through which a mixture containing the particles flows in 
order to have measurements performed thereon. 
The flow cell may be positioned to evaluate raw materials 
entering a process or to evaluate materials during a 

30 process . 



BRIEF DESCRIPTION OF THE DRJOTINGS 

The foregoing aspects and other features of the present 
invention are e3q>lained in the following description, 
taken in connection with the accompanying drawings, 
5 wherein : 

Fig, lA is a perspective view of a first form of an 
optical probe in accordance with the invention. 

Fig. IB illustrates a possible geometry of a view port 
10 and a probe for the optical probe of Fig. 1. 

Fig. IC illustrates focusing optics for the optical probe 
of Fig. 1 

15 Fig. 2 illustrates the change of diffuse reflectance 
response as an emulsion system breaks. 

Fig. 3 illustrates the change of diffuse reflectance 
response as a solid-particle system is milled. 

20 

Fig. 4 illustrates aui a^aratus in accordance with the 
invention for measurement of emulsion stability. 

Fig. 5 shows spectra formed by separation of sender and 
25 receiver fibers. 

Fig. 6 is a full 11-frame movie of sample subtracted from 
that of finished product. 

30 Fig. 7 is a fractional deviation movie. 
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Fig. 8 is an example of a frame where noise pushes the 
signal negative. 

Fig. 9 is a plot of fractional differences of sample C 
S movie and the mean movie. 

Fig. 10 is a plot of fractional differences for saxnple B 
as in Fig. 9. 

10 Fig. 11 is a plot of correlation functions for three 
creams 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Referring to Fig. lA, there is shown a perspective view 
15 of an apparatus 10 incorporating features of the present 
invention. Although the present invention will be 
described with reference to the various embodiments shown 
in the drawings, it should be understood that the present 
invention can be embodied in meuiy alternate forms of 
20 embodiments. In addition, any suitable size, shape or 
type of elements or materials could be used. 

Fig. lA shows the basic configuration of the system. A 
fiber optic or bimdle 2 directs light from a source 4 
onto a sample. There may be a transmission window or 
25 sight glass 6 that abuts the sample 8 to be studied. A 
second fiber optic or bundle 10 collects light from the 
sample and directs it to a spectrometer 12 and a detector 
14. In an equivalent configuration, the light entering 
through bundle 2 may already have been predispersed. The 
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spectrometer may be a photodiode array instrument such as 
those made by Ocean Optics or Micron Optical Systems, an 
AOTF such as those made by Brixurose, or an FTIR such as 
those made by PerkinElmer . 

5 

To minimize reflections from the glass 6, the glass may 
be treated with antireflection coating at the wavelengths 
used in the system. Alternatively, there may be a drop 
of index matching fluid 16 between the fibers and the 
10 glass. 

Another way to minimize back reflections is to rotate the 
angle at which light from the fibers meet the plane of 
the glass, as is shown in Fig. IB. 

15 

Varying the separation of the fibers will also tend to 
reduce spurious reflections at various distances between 
the fibers. 

20 

In yet another configuration, shown in Fig. IC, one or 
more of the fibers may be positioned with a focusing lens 
20 or a collecting lens 22 to allow imaging into portion 
24 or from portion 26 of the saxnple. The sainple may be 
25 in a flow cell 45. It is inqportant to note that 
variation of the focus position and/or the collection 
position of the lenses is another degree of freedom with 
which to examine the system, analogous to d. 

30 A sui^ort block 28 that holds the fibers is attached to a 
micropositioner 30, such as a single axis translation 
stage made by Newport Corporation, to move the fibers to 
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the desired separation distance d 32. The translation 
stage may be driven by a motor which can be interfaced to 
a computer to provide an integrated method of collecting 
spectral data at variable inter-fiber spacings, d. When 
5 d is close to zero, the system is in the backscattering 
mode, the standard diffuse reflectance configuration. 

Spectra Sd are collected by sets & with E consisting of at 
least one, but generally a plurality of spectra, at 

10 specified d values. Such sets are called ^^movies'' since, 
in general, they are snapshot spectra done in multiple 
positions or framers. For example, a set may consist of 
So, S2, S5, and Sio, where the subscripts denote the size 
of d in millimeters. Traditional diffuse reflectance 

15 measurements are a special case of E, containing only So* 



Data may be used in single*beam mode. In a preferred 
form, data may either be normalized by an absolute 

20 intensity as in the method of Schmitt and Kumar, or the 
data of E may be divided by So to observe deviations in 
the relative amplitude of the peaks of the various 
spectra. Thus, when the effect of increasing d is merely 
to uniformly increase the path length across the 

25 available spectrum, the new ratios, S^/Sq will all be flat 
lines . 

Since scattering depends on the ratio of the wavelength 
of the incident light X to the size R of the scattering 
30 particle (assumed spherical for simplicity) , when X/R » 
1, the particle is an isotropic scatterer. For large 
particles, i.e., when X/R » 1, the scattering tends to 
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beam in the forward direction. Two phenomena are 
anticipated: 

1. As X/R increases, more light is directed outward 
5 toward the large d values. This happens because the 
isotropically scattered beam has components going in 
every direction, including parallel to the interface. By 
contrast, large particles, which scatter in the forward 
direction, lack a significant component in the direction 
10 transverse to the direction of the injected light. Of 
course, it is the presence of multiple scattering that 
enables detection at significant d values in the first 
place . 

15 2. Within a spectrum, there is a chromatic effect. For 
example, for a bandwidth of 800 nm with X going from 800 
nm to 1600 nm there is a wide difference in X/R ratios. 
At the blue end, where X is smallest, particles appear to 
act larger than at the red end. These phenomenologies can 

20 be used to advantage in several ways: 



Particle Size and Coxuposition 

A series of E measurements may be made to represent the 
25 time evolution of a mixture as it is being mixed and or 
milled. In this process particles are frequently made 
smaller. Thus, the scattering properties of the mixture 
varies, or E is really E(t) where t is time. At some 
time tf the mixture will be completed. However, tf may 
30 vary due to a variety of external factors such as 
mixing/milling rate, initial size of particulates, etc. 
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If the E(tf) data is available, then a real time 
comparison of E(t) can be done. To obtain E(tf) a sample 
of the product that is considered to be of acceptable 
5 quality as determined by whatever standards are normally 
employed is separately measured. The quality standards 
will differ from product to product. For a food it might 
be texture, viscosity, or some other variable or set of 
variables. For a cosmetic it may be sheen or 

10 transparency. Nevertheless, each of these acceptably 
finished products can be linked with a E simply by 
measuring a sample. This E is by definition the E(tf) . 

In general, it will be possible to store a time evolution 
15 of E so that not only will it be possible to determine if 
a mixture is ready, but also how far along the process it 
is. 

Figures 2 and 3, which are representations of a S with 
20 only an So, i.e., a fixed spacing, backscattering probe, 
show that near IR radiation can be scattered back from 
particulate matter in liquid dispersions and that the 
resulting near IR spectrum provides information about 
both system caaposltlon and particle size. Additionally, 
25 the figures demonstrate that particles may be either 
self -forming droplets, as emulsions, or solid materials 
like minerals. 

30 It is universally known, that when water and an oil, like 
canola oil, are shaken together the two liquids form an 
emulsion. Basically globules of oil in water or water in 
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oil are formed. It is also well known, that the shaken 
system will eventually settle out and the oil phase will 
once again settle above the water phase. In the process, 
the moderate size emulsion droplets that had been created 
5 upon shaking, agglomerate. They then grow to larger 
droplets until the "droplet" size becomes the whole 
phase . 

The spectrum of the emulsion system chamges as the 
emulsion breaks, as seen in Fig. 2. These spectra were 
taken with an external Galileo probe using an indium 
antimonide detector. The Galileo probe is a fixed 
spacing, back-scattering fiber optic probe connected via 
the external port to a PerkinElmer FT spectrometer. 
Bright (100+) ink jet paper was used as a neutral 
reference background, so these reflectance spectra have 
units analogous to those for transmission spectra, i.e. 
percent reflected. The prominent water peaks around 1450 
nm and 1900 nm are clearly seen in the emulsion 
reflectance spectra. So, too are the structures of 
canola oil at 1700-1800 nm as well as in the 1200 nm and 
2300 nm regions. Thus, the system composition spectra 
can be recognized. 

25 But additionally, as a function of time, the spectra move 
in a downward direction. In other words, as the 
particles get bigger, the spectrum shifts downward. The 
number of particles is also decreasing since virtually 
all the oil in this 25% oil emulsion system is used up to 

30 make droplets. As the droplet size increases, the 
additional oil comes at the expense of sacrificing other 
particles. In this case there is also a moderate shape 
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change as evidenced by the difference in the spectra 
below 1350 nm and again between 1600 and 1900 nm. 

Samples of ink in the process of being milled were 
5 obtained by removal from the process wet-millers at 
different stages of COTipletion. Fig. 4 illustrates a 
downshift in the near IR spectnim in the course of this 
process . 

10 The ink goes from starting as a premix at the very top 
spectrum, through 2 manufacturing passes and then down to 
a final product as the lowest spectrum. 

Not all sample spectra move downward in time. Other 
15 data, not shown, exhibit a rise in spectral position as 
milling proceeds. 

The conclusion is that the near-IR spectrum, can in 
addition to showing details of chemical composition by 
20 the position of spectral peaks, also shows the relative 
particle size of mixtures of particles dispersed in a 
fluid. These particles could be liquid themselves, as 
with emulsions, or they could be solid materials as with 
the inks and dyes. 

25 

Emulsion Stability Analysis 

Emulsions are intrinsically unstable mixtures . It is 
frequently necessary to know the stability lifetimes for 
30 a variety of products. For excunple, the anticipated shelf 
life of creams and foods is in part determined by the 
phase stability of the product. Since the invention is a 
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sensitive detector of both scattering (relating to 
particle size) and absorption (relating to chemical 
composition) , it is useful in determining the status of 
an emulsion by moving the fiber optic probes to scan 
5 along the height hg 44 of a container with emulsion filled 
to height he 42, as shown in Figure 4. The light source, 
4, and the spectrometer, 12 are shown « 

10 The apparatus can be moved vertically along the column of 
a separating emulsion 40. Initially, all E will be 
virtually identical since the system is approximately 
uniform. As the system begins its phase separation, an 
interface will form, shown at height hi and E will vary. 

IS For example, if the intial mixture were a simple oil and 
water emulsion, the E data will evolve such that 
ultimately there will be a uniform set of them along the 
bottom section of the column, corresponding to the water 
phase, and a very different set of E values in the upper 

20 part of the column for the oil phase. However, 
considerably before full separation, it will be apparent 
that there is a discontinuous change in the E pattern. 
This discontinuity corresponds to the interface that is 
evolving . 

25 

In principle, only one measurement, taken at a constant 
spacing, d, which could be zero, may be needed. However, 
by taking multiple frames at differing inter-fiber 
separations, the analysis may be improved. 

30 

Average Particle Size And/Or Particle Size Distribution 
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Another objective of the invention is to determine 
average particle size and/or particle size distribution. 
Inasmuch as both the magnitude of light at increasing d 
values (item 1 above) and the relative amount of 
5 modification (or distortion) of the backscattered 
spectrum So as seen at higher d values (item 2 above) are 
a byproduct of particle size, there is data to be mined. 

Another feature to consider is that scattering behaves as 
10 X"'^ where n goes from approximately 4 in the Rayleigh 
(small particle) regime to < 1 in the Mie (large 
particle) regime . 

A way to process this data is to use samples of colloidal 
15 suspensions, preferably ones that can be diluted, and 
independently determine particle size. Next perform a 
multivariate analysis on E to extract the information. 
An example is a suspension of paint pigment, which unlike 
emulsion droplets, can be diluted in a routine manner. 
20 The diluted sample can then be measured with correlation 
techniqpies (dynamic light scattering) . 

Polydispersity or average particle size can be forward 
modeled from the equations in Scmitt and Kumar's work 

25 sinqply by summing up the results frOTi their calculations 
on monodisperse materials. Despite the fact that their 
calculation in equation 1 from the cited reference is the 
result of a random diffusion model, their results are 
dependent on entirely measurable quantities. Thus, in 

30 their modeling, they predict the reflectance of a system 
based on absorption coefficients, volume fractions of 
materials, surface reflectivities and, importantly, 

20 



scattering cross sections. These cross sections are 
themselves based on mean volume of spheres. Thus, for a 
measured reflectcuice and with all other properties but 
size known, a trial -and-error fit can be made to match 
5 the observed spectrum to the various trial particle 
sizes . 

Furthermore, the spectra calculated by Schmitt and Kumar 
are functions of the inter-fiber distances. By making a 
10 plurality of measurements at differing fiber spacings, a 
suite of spectra can be obtained. A best fit of the 
predicted spectra to the suite as a function of spacing 
using the estimated particle size would provide an 
optimal result. 

15 

This method could be extended to trial distributions of 
particles to fit a hypothesized particle size 
distribution to the measured reflectance. Ideally, this 
would be after selection of a size probability 
20 distribution and integrating over particle size. Thus, 
for example, if one assumed N different particle sizes, 
one could generate N diffusion equations and sum them up, 
suitably weighted, to give an overall prediction of 
spectral changes form polydisperse mixtures. 

25 

Utilizing a probe with a variable distance between the 
fdLber that injects light into the sample (the sender) and 
the fiber that collects the reflected light (the 
receiver) , results in different positions providing 
30 variable overall intensities. Fig. 5 shows a suite of 
such spectra for a finished commercial creeua. 
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When a fixed probe is used, comparison between different 
sanq>les can be achieved simply by overlaying their 
respective spectra and noting variations. The customary 
way to do this is to subtract one spectrum from the 
5 other. The closer the residuals are to zero, the more 
nearly alike are the two spectra, and therefore, the more 
nearly alike are the samples themselves. The analogous 
situation when comparing the suite of spectra, called the 
movie, is shown in Fig. 6. This, too, is a subtraction, 
10 but it generates a two-dimensional surface because the 
individual movie frames are plotted one behind the other 
in a 3-arreuigement. 

If the sample is finished product, its movie is the same 
15 as that of the movie of Fig. 6. The 3-D plot then 
collapses into a flat sheet at zero because all of the 
respective frames have matched spectra. In fact, the 
data at the higher frame numbers appears to be very close 
to zero, providing the impression that the two samples 
20 are well -matched at higher frame numbers . Such a 

conclusion would be misleading for the following reason: 
as Fig. 5 shows, the raw spectral data is much smaller at 
high frame numbers than at low ones. Consequently, 
subtracting two small numbers is bound to leave a small 
25 number as a remainder. For the early frames, where 
signals can exceed 30,000 counts, even a 3% variation is 
900 counts, a value larger than some of the spectra in 
the later frames. 

30 To improve the balance of the subtraction, a weighted 
difference is taken in place of the raw difference. 
In one implementation the two spectra in a particular 
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frame nuxniber are subtracted from each other, as before, 
but in addition the difference is divided by the spectrum 
of the finished product in a point-by-point manner. This 
scheme provides equal weighting to all the frames . 
5 Naturally, a more complex set of weights may be assigned 
if unequal weighting is appropriate. For example, a 
weighting that is proportional to the magnitude of the 
signal or to the square root of the signal might be 
chosen, where by ^'signal" the integrated area under the 
10 curve is meant. Other possibilities for defining signal 
could be considered. In effect, a fractional change in 
the spectrum is calculated. See Fig. 7. 

Under these conditions, the smoothness in the high-frame- 
15 number region is now replaced with a bumpy (positive and 
negative deviations) structure. The process of 

normalizing the data to provide a fractional difference 
rather than an absolute difference is seen therefore to 
improve the ability to discern variations in the data 
20 that are real and independent of frame number. 



This tool can be applied to a set of samples measured on 
the near IR instrument. The objective is to determine 
which samples appear to be most similar to other sauries 
25 and, correspondingly, which are most different. 

Since the objective is to find the degree of 
dissimilarity between samples, the best way to begin this 
effort is by taking maltlplB data sets on a slnglB 
30 saaple. Thus, three runs were performed on each 
particular sample, with the container moved to a fresh 
location between runs. Inasmuch as these data should be 
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very close to each other, this triplet sets a xainimum 
standard of variation that should be expected to be 
detectable . 

5 The contparison was xaade by: 

1. Taking the mean spectrum of all three runs for each 
frame . 

10 2 . Subtracting the mean spectrum from each of the three 
using the weighting previously described. 

3. Tabulating the results according to frame number. 

15 Table 1 shows reproducibility errors of 5 creams as a 
function of frame number. The table entries can be 
considered to be the overall percent error of mismatch 
between the specific run and the mean of the three runs. 
Thus, the errors are the difference between the mean of 

20 the three runs and the individual spectra. The larger 
deviation errors in the higher frames could be a useful 
indicator for unequal weighting of the subtraction data. 
For example, the weighting could be set to be inversely 
proportional to the deviation. 

25 



30 
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Table 1 



Sample Frame #: 6 7 8 9 10 11 Avg frame 10 

A 0.8665 1.6043 2.2022 4.0425 4.9299 21.5796 4.93 

0.9176 1.5954 1.8449 3.0223 5.1271 7.5346 

1.4925 2.4445 2.7322 4.3738 4.7347 4.8567 



B 



0.9038 2.1543 3.5576 26.1307 86.4859 311.7983 
0.8945 1.8154 18.7277 213.3659 254.6579 39.4301 
0.4202 1.3836 10.6072 13.2979 202.0071 16.6333 



181.07 



0.5055 0.6771 2.0096 5.4930 5.3456 30.5588 
0.4884 1.1762 1.5923 2.3089 4.3520 9.1865 
0.4840 0.7941 1.7240 4.6256 10.7439 7.7015 



6.81 



0.2636 0.4805 0.8893 1.1330 1.8702 8.2100 
0.2208 0.4094 0.9781 1.3273 3.0086 7.4759 
0.2577 0.6256 1.0410 2.7476 3.3683 5.5424 



2.75 



0.3231 0.3922 0.6993 2.6934 3.2948 48.9279 
0.2944 0.3910 0.5748 1.6568 1.5629 2.3534 
0.2388 0.4768 0.6467 1.4660 1.6531 34.5727 



2.17 



If frame 10 is (arbitrarily) chosen as the benchmark, 
samples A,C,D and E ajqpear to reproduce to about 7% or 
better, but senile B shows mismatches as high as 180%. 
This enormous error is fortuitous . The attempt to put 
all the frames on an equal footing by teJcing a ratio of 
the differences to one of the spectra has a built-in 
limitation at numbers near zero. Reference is made to 
Fig. 8 as eui example. 

An artifact is introduced in attempting to weight data 
points with very small numbers to the same degree as the 
larger euiplitude points. To reduce the effect of this 
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undesired situation, another technique is added: 
windowing. The region (window) from pixels 90 to 130 is 
cut out of the data cuialysis entirely. When this is 
done, the errors in all the points in table 1 go down in 
magnitude. The suxmoary for frame 10 can be seen in Table 
2. 

Table 2 

Full spectrum Windowed Spectrum 



A 4.93 4.26 

B 181.07 84.54 

C 6.81 1.89 

D 2.75 1.14 

E 2.17 0.48 



Three out of the five examples provide excellent 
repeatability when the bottom-scraping data points are 
eliminated. A fourth one, A, provides respectable 
results. B, though much improved is clearly an example 
of poor data. This suggests that the integrity of a 
sample data set should be inspected prior to moving on to 
another seunple. 

A closer look at Fig. 8 shows that at high bin values, 
(approximately 240) , there is another drop in the 
spectrum toward zero. A comparison between B and another 
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sample is illustrative of what a good quality control 
measure can provide. 

Fig. 9 is a plot of fractional differences of one movie 
(of three) and a mean movie. Referring in particular to 
5 saxnple C, most of the variation occurs in the region of 
the window excised, i.e. pixels 90 -130. The dynamic 
range is -4 to +8. 

In Fig. 10, one of the sample-B movies is generated under 
10 comparable conditions to those of Fig. 9. Not only is 
the amplitude of the variation in the 90-130 region 
dramatically higher for B, but there is also a huge 
variation in the 225+ region. No other samples have this 
problem. 

15 

It is possible that random instrument drift may account 
for some of the difference, However, variation in the 
second region was not accounted for, though it could have 
been. In other words, by adding more windows, problems 
20 of variability may be diminished. 

A reasonable general rule that would be helpful in 
ensuring that measured variations are the result of real 
spectral differences rather than fluctuations near zero 

25 is all regions that have data within an 8 of zero should 
be windowed out to the, data points that are n8. 
Naturally, c and n are particular to each system and are 
largely dependent on the noise. For example, a selection 
may be made to window out regions of the spectrum, which 

30 are less than 10 times the value of the noise. Clearly, 
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this could be a larger burden for the higher frames of 
the movie since, in general, their signal is lower than 
that of the early frames. 

Fractional differences are simply the difference of the 
individual run and the mean of the three runs divided by 
the mean of the three runs. 

Once the quality of individual samples is determined then 
variations between samples can be reliably q[uantified. 
The mean value of each of the samples is compared in 
Table 3. 

Table 3 





6 


7 


8 


9 


10 11 


AtoB 


1.8430 


2.8524 


5.2464 


9.5216 19.5637 


54.5614 


AtoC 


1.8642 


2.3544 


2.8646 


3.7874 4.3643 


4.8240 


AtoD 


1.2658 


1.7450 


2.2985 


3.8122 4.5601 


4.7147 


AtoE 


5.5110 


8.0512 


10.2515 


11.2953 12.5645 


13.6896 


BtoC 


0.6674 


1.1397 


3.3636 


7.7673 33.8393 


66.4533 


BtoD 


1.7980 


1.8870 


3.6698 


7.9811 26.8299 


193.8976 


BtoE 


6.9747 


10.1578 


13.7240 


17.0529 21.8182 27.6169 


CtoD 


2.0484 


1.8036 


1.8687 


2.4527 3.0835 


4.6691 


CtoE 


7.1682 


9.9882 


12.3888 


13.8554 15.1916 


15.8287 



DtoE 5.6766 9.0567 11.7461 13.7298 15.3005 16.3190 
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In accordance with the invention, one way to carry out 
reproducibility measurements during on-line operation is 
to allow the sample to flow. In this way, every movie is 
taking measurements on fresh san^le, not measured 
5 previously. This is the equivalent of moving the sample 
cup to examine different parts of the sample as done in 
the laboratory. Thus, a plurality of measurements, i.e. 
of movies, can be made during operation to ascertain the 
degree of reproducibility of in-line measur^ents. 

10 Naturally, this process is limited to those samples which 
themselves do not change appreciably during the 
measurement time of replicate movies. For example, if a 
movie takes 90 seconds to complete, three movies could 
take about 5 minutes • For a process that is completed in 

15 15 minutes, this duration is much too long. The sasiple 
would undoubtedly be altered during this time. In such a 
situation, single movies done either on the flowing 
sample, or preferentially on a static sample, would be a 
favored approach. On the other hand, for a multi-hour 

20 process, five minutes of sample observation is likely to 
be the equivalent of a single point in time, with a 
virtually constant coniposition . 



25 Mixture Homogeneity Detection 

When a diffuse reflectance spectrum is measured with 
sufficient S/N to constitute a repeatcd^le signal to 
within a specified tolerance, deviations from 
repeatability are indicators of a changing Seunple. The 
30 meaning of repeatable signal could be defined in at least 

29 



two ways. The simpler method would be to look at the 
overall intensity of the signal by integrating the area 
under the spectral curve. This is equivalent to having a 
single detector that does not break the scattered light 
5 into a spectrum, but looks at the intensity of the 
overall beam. A more sophisticated way would be to 
calculate the RMS difference of subsequent spectra to 
quantify the variation. This is a particularly useful 
method when the inhomogeneity derives from poorly 
10 dispersed constituents, generating spectral differences 
because of compositional variations. 

For an unchanging sample, there will be some variation of 
the diffuse reflectance spectrum attributable to random 
15 noise, instrument drift and other either controllable or 
uncontrollable variables. It is vital to know the 
repeatability limits of the instrument. 

When measurements are made showing non-repea table 
20 results, i.e. spectra outside the expected limits of 
error of the instrument, these results are indicative of 
a change in the sample. Since most on-line measurements 
of mixing/milling processes are done on flowing systems, 
the spectral result is really due to an average over a 
25 volume of material. For exaznple, a near-IR fiber optic 
probe could be on a portal in a vat looking into a 
stirred stream; directly immersed into the stirred 
stream; or, at a sample cell that is part of a slip 
streeun whereby some of the dispersion from the vat is 
30 pumped out and through the cell and generally allowed to 
return to the vat. 
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During the measurement interval a volume proportional to 
the measurement time passes through the detection region. 
As the detection time is increased, more fluid flows by 
and greater amounts of the dispersion are averaged. 
5 Conversely, shortening the measurement time is equivalent 
to looking at smaller volximes of the dispersion. 

Now it is typically true that the larger the volume 
examined, the greater the degree of apparent homogeneity. 

10 Items present in low abundance, particularly large 
aggregates will appear to fluctuate greatly in number 
when small volumes are chosen. This result is a simple 
consequence of the standard deviation being proportional 
to the square root of the number of elements, N, and the 

15 relative deviation is just the inverse of this quantity 
or N~^^^. As N gets larger, the relative error decreases. 

A common way to manually assess the homogeneity of a 
dispersion is to remove a sample and place it between to 
20 microscope slides. Then it can be held to a light and 
examined for uniformity. 

An improved, automated way to achieve the same goal, is 
to examine the near-IR (or indeed any other radiation 

25 that can be scattered from the dispersion) as a function 
of integration time of the measurement. For example, if 
a measurement is done over a 15 second interval, sampling 
the voliame that traverses the probe in that time, the 
same type of measurement may be done over briefer 

30 intervals, say 10 seconds, 5 seconds, 1 second. The only 
criterion is that the repeatability of the system at 
these different integration times must be kept in mind. 



These criteria had previously been established with 
unchanging sauries. 

Thus, the homogeneity of the system can be derived from 
5 the time scale of integration, which is equivalently 
proportional to the volume of material, needed to obtain 
consistent results within the criteria of 
reproducibility . 

10 The mixing/milling process frequently has multiple 
requirements including full dispersion of the chemicals 
and a grinding down of particle size. The time scales 
for each of the activities must be kept in mind. For 
example, the system may have rapidly become homogeneous 

IS but is undergoing constant and significant particle size 
reduction. This event, however, would be detectable by 
the consistent (monotonic) change in the series of 
spectra. An inhomogeneous system would have, on average, 
as many upward and downward fluctuations in the spectral 

20 data. 

The quality of hranogeneity of a dispersion undergoing a 
particular manufacturing process can be quantified by 
looking at the time scale needed to achieve repeatability 

25 of the signal. When doing so, it is iniportant to make 
sure that flow rates are taken into consideration since a 
faster slipstream for example, would mean a greater 
amount of material per unit time and would give the 
illusion of greater homogeneity. It is good practice to 

30 fix the flow conditions to a specified reference standard 
before making conclusions about uniformity of the 
dispersion. 
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For a manufacturing process done according to a fixed 
protocol, the integration time procedure is a useful 
method of judging how close a product' s uniformity is to 
5 that of previous successful products. 

Validation of Coinponents 

A second measurement port may alternatively be used in 
the mixer either instead of or in addition the probe that 

10 observes the mixture. This would be a dual purpose cell 
capable of measuring either in transmission or 
reflectance to examine incoming material and compare its 
spectrum to a database of known materials. In one 
embodiment, the sample cell would be multiplexed to make 

15 either type of measurement. In the other, only a 
transmission or a reflectance measurement would be 
needed. 

For example, in the manufacture of dispersions comprised 
20 entirely of transparent liquid materials, only a 
transmission cell would be needed. In the case of a more 
conqplex mixture involving powders, a dual cell would be 
appropriate, or just a diffuse reflectance cell for the 
case that only the powders are of concern. 

25 

The cell may have a single light source and a fiber or 
other optical method to couple diffuse reflectcmce from 
the front surface. It would also be equii^ed with a 
fiber at the rear of the cell opposite the illumination, 
30 for transmission measurements. In another enibodiment, no 
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fiber optics would be necessary; a detector would be 
inline . 

Another way to inclement this cell would be in series, 
5 where the first cell is in transmission mode and the 
second is in reflectance, or vice versa. 

As material is added, the spectra are taken. For a given 
10 flow rate measured with an auxiliary device, the amount 
of material inserted can be determined by the duration of 
time that the spectnmi of the particular material is 
measured. Alternatively, a weighing device may be used. 

IS A ccMoputer or similar storage device records this 
information for later retrieval to validate the addition 
of proper cc»rponents. 

In a preferred embodiment, the full mixture is monitored 
20 with the near-IR, as well, to verify tiiat that the 
particle size and overall conrposition of the mixture is 
as required. 

Laser Method to Exaaiine Creams 

25 

Fig. 11 shows autocorrelation data from laser light 
scattering in three cream samples . Notwithstanding that 
the classical diffusion theory is suspended for 
interpretation of the data, the fact that reproducible 
30 correlation functions are obtainable is significant. 
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The ^jective is to be able to associate particular 
correlation functions with specific chemical 
compositions. In this way, the correlation-data strategy 
resembles much of the earlier discussion of broadband 
5 analyses. Both methods rely on the prior knowledge of a 
defined goal material, and both accumulate standard 
results which can be stored and used later for ccnnparison 
with new saxoples. This association of a particular shape 
correlation function with a specific product or with a 

10 product in a stage of being processed is useful for 
making distinctions among products with similar 
appearance. More importantly, it can be used as a set of 
curves representing a cream or lotion during 
manufacturing. The correlation measurements can be made 

15 on the samples, preferably on-line, during mixing/milling 
to determine the current state of the system. By 
comparing the correlation functions to a known set taken 
as a library reference, the operator can judge the 
current condition of the process. The process can then 

20 be terminated when the curve corresponding to the 
finished saxiple is obtained. 



The raw correlation function is subject to variations 
25 that significantly modify its appearance, but in ways 
unrelated to the chemical system. These artifacts need to 
be eliminated. First, subtract the theoretical baseline 
from all of the samples. Correlation functions have 
baselines equal to the product of the total number of 
30 counts times the square of the average number of counts 
per channel. Then consider the possible statistical 
variation in the laser coherence, a phenomenon that is 
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likely to vary significantly in very short time 
intervals, maybe even as short as a second. This 
fluctuation impacts the overall amplitude after baseline 
subtraction. When this fluctuation is taken into 
consideration simply by normalizing the remaining 
correlation function, the results shown in Fig. 11 are 
obtained. This figure shows data taken in triplicate for 
cream samples in three states of manufacture. 

There are numerous ways to characterize the evolving 
correlation functions. One way is to use the standard 
cumulant analysis, which after the filtering already done 
is basically a polynomial description of the curve. This 
can be done to arbitrary order, though in practice it is 
rarely carried out beyond fourth order. 

A second way would be to assemble a series of samples and 
measure their correlation functions and then perform a 
cluster analysis so that particular formulations would be 
associated with a given cluster. A measured unknown 
could be assessed for membership in or proximity to a 
given cluster. 

It should be understood that the foregoing description is 
only illustrative of the invention. Various alternatives 
and modifications can be devised by those skilled in the 
art without departing from the invention. Accordingly, 
the present invention is intended to embrace all such 
alternatives, modifications and variances, which fall 
within the scope of the appended claims. 
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